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The centrality dependence of elliptic flow and how it is related to the 

physics of expansion of the system created in high energy nuclear collisions is 

discussed. Since in the hydro limit the centrality dependence of elliptic flow 



D . is mostly defined by the elliptic anisotropy of the overlapping region of the 



colliding nuclei, and in the low density limit by the product of the elliptic 

anisotropy and the multiplicity, we argue that the centrality dependence of 

^ . elliptic flow should be a good indicator of the degree of equilibration reached 

f— ^ , in the reaction. Then we analyze experimental data obtained at AGS and SPS 

^D . energies. The observed difference in the centrality dependence of elliptic flow 

- could imply a transition from a hadronic to a partonic nature of the system 

evolution. Finally we exploit the multiplicity dependence of elliptic flow to 

^ , make qualitative predictions for RHIC and LHC. 



H ' The goal of the ultrarelativistic nuclear collision program is the creation of the QGP - 

quark-gluon plasma - the state of deconfined quarks and gluons. It is understood that such 
a state requires (local) thermalization of the system brought about by many rescatterings 
per particle during the system evolution. It is not clear when and if such a dynamical 
thermalization can really occur. An understanding of these phenomena can be achieved 
by considering elliptic flow [|l| recently studied at AGS and SPS energies. It will 
be shown how the centrality dependence of the strength of elliptic flow, V2, deflned as the 
second coefficient in the Fourier decomposition of the particle azimuthal distribution IQ, is 
an indicator of the degree of equilibration (thermalization) achieved in the system. 



Our qualitative conclusions are based on the observation, that in the hydro limit (which 
we equate in our discussion to complete thermalization) and in the opposite limiting case, 
the low density limit, (where dynamical thermalization is not expected) the centrality de- 
pendence of elliptic flow is different. In the hydro limit, the mean free path is much less than 
the geometrical size of the system. The centrality dependence of flow is totally governed in 
this case by the initial geometry (eccentricity), the latter being roughly proportional to the 
impact parameter. In the low density limit, the mean free path is comparable to or larger 
than the system size. The final anisotropy in this case should be proportional to the ratio of 
the system size to the mean free path (the number of collision) . The anisotropy vanishes in 
the limit of infinite mean free path. The latter in its turn depends on the particle density, 
which is largest for central collisions and vanishes for very peripheral collisions. Note that 
the factors involved change drastically with centrality. One could imagine other reasons 
for centrality dependence of elliptic flow in the hydro model, such as the initial conditions, 
viscous corrections, resonances, or effective volume corrections, but we expect that all these 
other factors have a much weaker dependence on the impact parameter. By considering 
the two limiting cases we hope to highlight qualitative considerations important for under- 
standing the degree of thermalization and the partonic or hadronic nature of the collisions. 
In essence we present a framework for examining these questions experimentally, however 
at the moment, it is mainly the experimental data which are not adequate to answer these 
questions convincingly. 

We will often use the term "physics of the collision" . By this we mean both the degree 
of equilibration and whether the hadronic picture in terms of nucleons, pions, etc., or the 
partonic picture in terms of deconfined quarks and gluons, is more applicable to the evolution 
of the system. The partonic picture in our view is similar to a QGP but the system is not 
necessarily thermalized. 

Low Density Limit 

To discuss the centrality dependence of f 2 more quantitatively, we start from the hypoth- 
esis that the system is not dense and its evolution can be described by the first correction 



to the collisionless limit 0. Physically this means that the rescattering occurring during 
the system evolution changes the particle momenta very little on the average and the corre- 
sponding change in the distribution functions can be treated in first order as perturbations. 
Under this assumption the final elliptic flow, f2, is proportional to the initial overlapping re- 
gion elliptic anisotropy, e, (introduced in flow analyses in []I| and in its present form in [^,0) 
and to the initial particle space density which defines the probability of particles to rescatter 

i- 

The initial geometry of the overlapping zone can be evaluated in a simple Glauber type 
model with a Woods-Saxon nuclear density. The results are weakly dependent on the weights 
used [0]. What is important is that if one wants to compare different energies, e.g. AGS, 
SPS and RHIC, the nuclear geometry cancels out, and only the dependence on multiplicity 
is left. This is true provided that the "physics" of of the system evolution stays the same. If 
it changes then the scaling with multiplicity will be violated. This is a very important point 
if one reads it the other way around: if scaling is not observed then probably the physics 
has changed. 

Under the assumption that the system is relatively dilute the momentum anisotropy is 
proportional to the spatial anisotropy, but also the particles must scatter to probe that 
anisotropy. Thus, the spectra distortion is directly proportional to the spatial anisotropy 
and the number of rescatterings, or the particle density in the transverse plane. In this limit 
the final elliptic flow (see a more detailed formula in [^) 

IdN 

where S = nRxRy is the area of the overlapping zone, with i?^ = (x^) and B^ = {y"^) 
describing the initial geometrical sizes of the system in x and y directions, respectively. 
(The x-z axes lie in the reaction plane). The averages include a weighting with the number 
of collisions along the beam axis. The initial space elliptic anisotropy is defined as 

y ~ ^ / \ 

it^ -\- Thy 



In our calculation we use a Woods-Saxon parameterization of the nuclear density with pa- 
rameters Ra = 1.12 ■ y4^/^, and a = 0.547 fm. More information on the effect of different 
weights and the values of R^, Ry, S and £ as a function of impact parameter can be found 
in [^. The proportionality coefficient in Eq. (1) is defined by the "physics" of the rescat- 
tering. If the physics is the same in central and peripheral collisions then Eq. (1) yields the 
centrality dependence of f2. 

Hydro Limit 

As follows from Eq. (1) the elliptic fiow increases with the particle density. Eventually 
it will saturate P] at the hydro limit, which would mean complete thermalization of the 
system. In this regime the centrality dependence of elliptic fiow is mainly determined by the 
initial elliptic anisotropy of the overlapping zone in the transverse plane , and the ratio 
of the two should be approximately constant as shown in the first such calculations done 
by Ollitrault [|l|. From his results it follows that {v2/e)hydro ~ 0.27 — 0.35, depending on 
the equation of state used (with or without QGP)[[ The calculations give a somewhat 
smaller fiow, resulting in {v2/e)hydro ~ 0.21 — 0.23 (partly due to the realistic treatment of 
resonances which decrease the pion fiow by about 15%). Note that in both calculations, |T| 
and [^], the longitudinal expansion of the system is treated analytically assuming Bjorken 
scaling. Real 3D hydro calculations would be very useful, although we do not expect that 



^To avoid confusion, note the difference in definitions of e used in Eq. (2) of this paper and a^ 
from ||l|. For Pb+Pb collisions the maximal value of e ~ 0.44 compared to a ~ 0.3. . Then, the 
results [y yield V2 /e ~ 0.55 — 0.7, where V2 means the elliptic flow weighted with pf. Recent 

f 2 "I 

calculations Q show that the particle elliptic flow is related to this quantity as V2 ^ O.bvl ■ 



they would greatly change the centrality dependence. 
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FIG. 1. Top: comparison of elliptic flow, V2, for pions from RQMD ver. 2.3 (filled circles) with 
the dependence expected for the low-density limit (solid line) and that expected for the hydro limit 
(dashed line). Bottom: ratios of W2/v2'^^5 and V2/v2^^^^ . 

RQMD 

Before discussing the experimental data we will first consider a realistic model. We 



take RQMD v2.3 |T0[ for our calculations. Fig. 1 top shows the comparison of the directly 
calculated V2 of pions in Pb+Pb collisions at 158 GeV-A collisions at mid-rapidity (—1 < 
y < 1) with the expectation from the low density limit, V2^^ (Eq. (1) normalized to the 
same area under the curve in order to illustrate just the centrality dependence.) One can 



see rather good agreement, which suggests that RQMD is close to the low density limit 
even as one scans the centrality from peripheral to central collisions. (In this version of 
RQMD no QGP is simulated.) This is not that striking a conclusion, considering that no 
hydro-type behavior has ever been observed in RQMD. Note that the low density limit does 
not mean a low number of total rescatterings. The number of rescatterings can be large 
provided all of them are relatively soft and the particle momentum changes little compared 
to the initial momentum. The cross section which enters the equations is the transport 
(not total) cross section (see f^). The centrality dependence expected for the hydro limit is 
shown on the same plot by a dashed line also normalized to the same area under the curve 
{v2^^^'~' ~ 0.0595). Note the large difference between the two curves, which was not noted 
in 1^. Fig. 1 bottom shows that the ratio of V2 to the expected functional form is flat for 
the low density limit but not for the hydro limit. A centrality dependence similar to the low 
density limit was also observed in [^ where a computer simulation of a pion gas expansion 



was studied. 
Data 

Now let us turn to the experimental data. At AGS energies the elliptic flow of charged 
particles and of transverse energy was measured by the E877 Collaboration. Unfortunately, 
the publication P] containing the detailed pseudorapidity dependence for each centrality 
lacks a figure showing just the centrality dependence. Our estimates based on their data p| 



of charged particle flow at midrapidity are presented in Fig. |[ 
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FIG. 2. Elliptic flow at the AGS (open circles) and the SPS (filled squares). 



The data indicate that at AGS the flow peaks at niid-centrality|, consistent with the 
low density limit prediction and no change in physics with centrality. At this energy some 
decrease of elliptic flow in peripheral collisions can be also attributed to shadowing by 
spectator matter. At SPS [|13|, preliminary data indicate that the elliptic flow peak moves 
towards peripheral collisions. This fact itself would hint at the hydro-dynamical picture of 
the system evolution. A more detailed look at the data shows that this is unlikely. First, 
the maximal value of elliptic flow {v2 < 0.04) is signiflcantly less than predicted by hydro 
calculations |[l|,^ (about 0.09-0. 1)|. Second, in the hydro limit elliptic flow should depend 



^A similar centrality dependence of transverse energy flow (from the same data M) can be found 
in the thesis of Chang [|12| . 

^In Q agreement was claimed between hydro and the NA49 mid-central data |^ leading to their 
conclusion of complete equilibration. However, this comparison was done for pt < 0.3 GeV/c and 



only on the initial space elliptic anisotropy, e. The preliminary NA49 data indicate that the 
ratio V2/S, at least for semi-central collisions, is likely increasing with centrality [|13| (see 
the data presented in Fig. 3 below). This centrality dependence (natural for the low density 
limit) implies that we still could be far from the hydro regimeQ. 

Assuming that at SPS the hydro regime is not reached yet, the observed centrality 
dependence of elliptic flow would indicate that the physics of the system evolution is different 
in central and peripheral collisions. Elliptic flow peaks at more peripheral collisions because 
the central collisions exhibit too little flow compared to that expected from the AGS data 
scaled with multiplicity. A natural explanation for this would be that peripheral collisions 
are described by hadronic (re) scatterings (the same as at the AGS in both peripheral and 
central collisions) while in central collisions partonic physics becomes important. One of 
the possible mechanisms responsible for the change could be a color percolation occurring 
at high parton densities in the central collisions and discussed in more detail below. 

Discussion 

Summarizing, our view of the overall picture is: at AGS energies, the physics of rescat- 
tering which defines the system evolution is hadronic in nature, while at SPS it is the same 
for peripheral collisions, but for central collisions the physics is likely to be partonic. The 
partonic picture will remain at RHIC energies, with some extension toward more peripheral 



it could be that the pt dependence of V2 in the hydro model does not agree with experiment. 
^ One can argue that, taking into account systematic uncertainties, the preliminary SPS data for 
V2/£ are consistent with being constant as a function of centrality. In this case it would indeed 
mean that the system has equilibrated and the hydro regime has been reached. The low absolute 
strength of the elliptic flow in this case would indicate that the equilibration happens at a rather 
late time when the spatial anisotropy e has decreased due to initial "free streaming". We do not 
exclude this possibility but must wait for the final SPS data and the coming RHIC data to answer 
the question. 



collisions. At RHIC equilibration becomes more important, but it is not clear if complete 
thermalization will be reached. At LHC energies the parton densities could become so high 
that (partonic) rescattering would lead to dynamical equilibration of the (partonic) system 
(creation of regions of real QGP) and consequently to a hydro-dynamical type of system 
evolution. 

The above picture for collisions of heavy nuclei implies that the shape of the centrality 
dependence of elliptic flow would change continuously with beam energy. At AGS, the elliptic 
flow is peaked at an impact parameter value slightly higher than Ra-, just as prescribed by 
the low density limit. At SPS energies the peak moves toward more peripheral collisions 
because possibly the physics of relatively central collisions may have changed from hadronic 
to partonic, which leads to weaker flow than one would expect taking into account the 
increased multiplicity. If thermalization is not reached at RHIC, the elliptic flow peak 
could move back toward mid-central collisions because the physics of the peripheral and 
central collisions will be the same - partonic rescattering, unlike the situation at SPS when 
peripheral collisions are driven by the hadronic rescatterings resulting in relatively large flow 
signal. At even higher energies at LHC, the elliptic flow should peak at more peripheral 
collisions just as predicted by hydrodynamic calculations. 

The schematic overall picture based on these observations is presented in Fig. ^, where 
the ratio of elliptic flow to the initial space elliptic anisotropy is presented as a function of 
initial particle density.^ In this plot we use the experimental charged particle multiplicity, as- 



^At the moment this plot is qualitative as many things shown have large uncertainties. The hydro 
limits can depend slightly on the initial particle density [||,P] and, more importantly, on the time 
of thermalization of the system. The values shown are an average of the results of HJ^. The 
predictions for the case without QGP are only for the EoS of a massless pion gas. Resonances 
can soften the EoS and lead to weaker flow. The uncertainty in the experimental points is mainly 
from the determination of the collision centrality required for calculation of the initial space elliptic 



suming that it is proportional to the total particle multiplicity and also to the initial particle 
multiphcity. For the experimental values we use dNch/dy at mid-rapidity from n,|T5|. 
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FIG. 3. Elliptic flow divided by the initial space elliptic anisotropy at the AGS (open circles) 
and the SPS (filled squares). The shaded area shows the uncertainty in the SPS experimental data 
due to the uncertainty in the centrality determination. See text and footnote for the description 
of the curves and hydro limits. 

In the limit of very low density the objects which rescatter must be hadrons. At some 
critical density a partial deconfinement happens. Parton density becomes high enough such 
that the color parton can propagate in the perpendicular plane without hadronization. Each 
parton is always close enough to other partons which screen its colorQ. Once the motion in 



anisotropy and the area of the overlapping region. The data points correspond to the centrality 
determined from the fraction of the total cross section corresponding to each centrality bin. Higher 
centralities were estimated from experimental measurement of the number of participants [0|. Fi- 
nally, the smooth dashed curves are just schematic illustrations for hadronic and partonic scenarios 
and the solid curve includes a transition between the two. 
^This picture is very close to the deconfinement (color percolation) model discussed by Satz |1G] 
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the perpendicular plane becomes easier (there is no need for hadronization) , the elliptic flow 
decreases. Note that the system still can be far from being dynamically thermalized, which 
would occur only at even higher particle densities. Even more important, such a significant 
change in the behavior of V2/S can only happen if the system is not thermalized. See also the 
discussion of this question in [^,^ along with the discussion of the possibility of observation 
of the QGP to hadron gas phase transition. 

To prove or disapprove the picture described above one needs more accurate data on the 
centrality dependence of elliptic fiow. We would like to emphasize the importance of flow 
measurements not only at medium impact parameters but in the full range of centrality 
including rather central collisions where the anisotropic flow is small. The measurement of 
elliptic flow and its centrality dependence at RHIC thus becomes very important. Different 
models predict different rapidity densities for RHIC and LHC Assuming that they are higher 
than at SPS by factors of 2 and 8, respectively, we have indicated the regions expected for 
Au+Au (Pb+Pb) collisions in Fig. ^. The measurements of elliptic flow in collisions of 
lighter systems (e.g. Cu+Cu) are also very important since they would cover the region of 
the SPS Pb+Pb data and would be useful in testing the above picture. The new SPS data 
taken at 40 GeV-A energy are also of great interest since they would bridge the two other 
sets of data and may scan the onset of deconflnement from hadronic to partonic physics. 

Note that our picture of nuclear collisions and QGP production is different from what 
is usually discussed, which assumes thermal equilibrium even at rather low beam energies, 
when QGP is not expected, and then with an increase in collision energy, formation of 
regions of QGP. We believe that what could happen is that the deconflnement can occur 
before dynamical thermalization is achieved Jl^ and that the centrality dependence of elliptic 



flow would be a good indicator of this. 
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